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ABSTRACT 

We consider a synchronized, circular-orbit binary consisting of a polytrope with in- 
dex n and a point-mass object, and use a self-consistent field method to construct the 
equilibrium structure of the polytrope under rotational and tidal perturbations. Our 
self-consistent field method is distinct from others in that the equilibrium orbital angu- 
lar velocity is calculated automatically rather than being prescribed, which is crucial for 
obtaining apsidal motion rates accurately. We find that the centrifugal and tidal forces 
make perturbed stars more centrally condensed and larger in size. For n = 1.5 polytopes 
with fixed entropy, the enhancement factor in stellar radii is about 23% and 4 — 8% for 
fj- = I and ~ 0.1 — 0.9, respectively, where n is the fractional mass of the polytrope 
relative to the total. The centrifugal force dominates the tidal force in determining 
the equilibrium structure provided fi ^ 0.13 — 0.14 for n ^ 1.5. The shape and size of 
rotationally- and tidally-perturbed polytropes are well described by the corresponding 
Roche models as long as n ^ 2. The apsidal motion rates calculated for circular-orbit 
binaries under the equilibrium tide condition agree well with the predictions of the 
classical formula only when the rotational and tidal perturbations are weak. When 
the perturbations are strong as in critical configurations, the classical theory underesti- 
mates the real apsidal motion rates by as much as 50% for n = 1.5 polytropes, although 
the discrepancy becomes smaller as n increases. For practical uses, we provide fitting 
formulae for the density concentration, volume radius, coefficient of the mass-radius re- 
lation, moment of inertia, spin angular momentum, critical rotation parameter, effective 
internal structure constant, etc., as functions of /x and the perturbation parameters. 

Subject headings: binaries: close - method: numerical - stars: interiors - stars: rotation 



1. Introduction 

Binary systems are an ideal laboratory for the study of stellar astrophysics. They not only 
provide direct information on stellar mass, radius, and luminosity, which are all fundamental for the 
theory of stellar evolution, but also allow to probe stellar structures by measuring apsidal motions 
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affect the rotation rate and internal structure of the other component, one expects that the stellar 
quantities such as radius and effective ter nperature of b i nary s tars inay dif fer from those of single 



stars, as confirmed by observations (e.g., iMalkovl l2003l . 120071 : iRibasI 120061 ). While the literature 



abounds with studies on the rotational and tidal distortions of stars, most of the studies are either 
approximate or based on prescribed rotation profiles, as detailed below. In this work, we employ a 
fully self-consistent method to investigate the quantitative effects of centrifugal and tidal forces on 
equilibrium polytropes in synchronously-rotating, circular-orbit binaries. 

Finding e quilibrium confi g uration s of polytropic stars under rotational (and/or tidal) distortion 
dates back to Chandrasekhar ( 1933a b 3) who used a perturbation theory in which an equilibrium 
density distribution was regarded as a small departure from a corresponding unperturbed Lane- 
Emden function. Chandrasekhar considered only first-order terms in the rotation parameter v 
(see eq. [17] for definition), corresponding t o slow rotation. Later, this perturba tion method was 



extended to include terms of second order (jAnand 



197 



1981 



1968 



Geroyannis et al.lll979l ) and thir d order 



Gerqvannis &: Valvil Il984l ). The first-order perturbation theory was further refined by iLinnell 



who retained the first-order terms in the potential expansion, while allowing for large 
deformation of th e equilibrium configurations from the unperturbed polytropes. More recently, 
Gerovannisl (|l988l ) showed that the first-order perturbation theory combined with "complex Lane- 



Emden functions" beyond their first roots gives much improved results for rotating polytropes. 



Kippenhahn &: Thomas I (jl970l ) used a quasi-spherical method in which (1) the rotation ve- 
locity is assumed to be spherically symmetric rather than being cylindrically symmetric and (2) 
only the radial component of the centrifugal force is consi dered in the force balance; these assump- 
tions make rotating stars essentially spherical. Recently, lArbutinal (|2009l ) used a similar method 
to study the change of the total angular momentum due to tidal perturbations and the associated 
tidal instability. While the quasi-spherical method was able to handle the effects of rotation on 
stellar evolution in a very simplified way, it readily fails when the stellar shape becomes quite 
non-spherical. 



To study the internal structure of binary stars, iKopall (j 19721 ) constructed the Roche model 



where the actual equipotential surfaces of disto rted stars are approximated by the corresponding 
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with very high central concentration, it assumes Keplerian orbits and thus cannot give any clue to 
the apsidal motions of binary stars. 



On the other hand. lMonaghan &: Roxburghl (j 19651 ) developed a so-called "double-approximation 
method" in which a distorted star is divided into two parts: a core containing most of the stellar 
material and a non-self-gravitating, low-density envelope. While the envelope is highly susceptible 
to applied centrifugal and tidal forces, the distortion of the core is assumed to be small even in 
the case of critical rotation. Depending on the method to calculate the distortion of the core. 
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the double-app r oximation method has seve ral di fferent branches: iMonaghan &: Roxburghl (|l965l ) 



i2pi 

JacksonI (|l97d ). iDurnev &: Roxburghl (119701) . andlNavlor fc Anand 



(Il972allbl) use d first- order per 



turbation theory of IChandrasekhaij (Il933all: iMartinI (|l97d) and ISiiigh &: Singhl (|1984| ) included 
terms up to of second order; Mevnet Sz Maedei (1997) and Lai et al. ( 20od ) employed the quasi- 
spherical method. 

Although the aforementioned work have improved our understanding of the rotational and tidal 
effects on stellar structure, these are o nly approxirn ate and do not yield accurate results when the 
perturbations are rather strong. It was I JamesI (|1964| ) who first sought for direct numerical solutions 
for rotationally distorted polytropes without any assumption. By expanding the potential-density 
pairs in Legendre polynomials up to of tenth orde r and integrating a set of the resulting ordinary 
differential equations simultaneously, I JamesI (119641 ) obt ained solutions, w i th tru ncation errors less 
than 0.2%, for rotating polytropes. On the other hand. lOstriker &: MarkI (jl968l ) introduced a self- 
consistent field (SCF) method that provides a concise, accurate, and very efficient computation 
scheme. In the SCF method, the equilibrium density distribution is found by solving the Poisson 
equation as well as the equilibrium condition, alternatively and iteratively, for the potential-density 
pairs. The SCF method is so powerful that it has been apphed to various problems including 
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the st ru cture of rotatiii g stars (jOstriker &: Bodenheimed Il968l : iMarkI 1 19681: 



1986al) and galactic dynamics of coll isionless svst ems (iHernguist fc Ostrikei 



Weinberg 



199S 



Clement 



1974 



Stabler 



Gnedenll2003l ). Recently, iMacGregor et al.l (|2007l ) used a 



1992 



modified version of the SCF method to construct stellar models fo r differentially-rotat ing, low- 
mass, main-sequence stars with a realistic equation of state (see also I Jackson et al.ll2005l ). 



Wh il e most SCF methods for distorted stars are two-dimensional under rotational symmetry. 



Hachisul (|l986bl ) formulated a three-dimensional SCF method to calculate structures of rotating 
multi-body systems in which tidal effects are implicitly included. He showed numerically that 
there is an equilibrium sequence along which an ellipsoid takes on a dumbbell-like shape as the 
angular r nomen tum of the whole system increases, eventually separating into a binary. Although 
Hachisul (jl986bl ) calculated the density structures and shapes of multi-body systems, these are for 
equal-mass objects (owing to the impos ed azimuthal symmetr y) with a few prescribed laws for 



spin and orbital angular velocities. Later. iHachisu et al.l (Il986al |bl) generalized the method to allow 



for unequal-mass binaries, but the rotational velocity is still given externally rather than being 
calculated self-consistently, which makes it impossible to obtain information on apsidal motions. 

Finding accurate internal structures of distorted stars in close binary systems are also of great 
importance in testing the classical theory of apsidal motions. Apsidal motions arise when the 
gravitational force due to a star is different from that of a Newtonian point-mass. The chief cause 
of perturbations to Keplerian orbits is the rotational and tidal distortions that make the figures of 
binary stars deviate from spherical symmetry. Also, the separate or combined effects of the general 
relativistic correction, the pres ence of a third body, and non-equilibrium tides can be al so non- 
neglig ible (see recent review of ICimenea l2007l ) . For binary stars in near-circular motions, iRussell 
(jl928l ) showed that the advance rate, vj, of the apsidal line due to the tidal and rotational distortion 
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of star I with mass Mi caused by the attraction of star II with mass M2 is given by 

5 




(1) 



where i^x = {G[Mi + M2]/A^y^'^ is the Keplerian angular velocity, i?i is the mean radius of star 
I, A is the orbital separation, and k2 is the dimensionless apsid al mot i on con stant. Through more 
rigorous, first-order perturbation analyses, ICowlind (|l938l ) and ISternI (119391 ) generalized equation 
([T|) to make it applicable to binaries with eccentric orbits, which has been widely used to estimate 
the degree of the internal density concentration via k2 from observations. 



However, the classical formulae of iRusselll (119281 ) . ICowlingj (jl938l ). and ISternI (jl939l ) have two 
limitations. First, by ignoring the cross terms of spherical harmonics in expanding equipotential 
surfaces, they are valid in a strict sense only when the tidal and rotational distortions are weak. 
Second, by assuming an instantaneous adjustment of a star to time-dependent tidal force, they fail 
when the orbit al period is shorter th an t he oscillation periods of tid ally excited modes inside the 
perturbed star. iQuataert et al.l (jl996l ) and lSmevers Sz Willemd (j200ll ) relaxed the second limitation 
to study the effects of dynamical and resonance tides (still based on a linear perturbation theory), 
finding that the classical formula e are accurate as lorig as t he ratio of the modal to orbital frequencies 
is less than about 10 (see also IClaret &: Willema |2003| ) . The validity of the classical formulae 
is nonetheless questionable for close binaries, especially in the critical configurations, where the 
distortion of a stellar shape is large enough for nonlinear terms to be non-negligible. 

Despite a long history of research, therefore, studies are still lacking as to how the stellar struc- 
ture, shape, size, etc. of a distorted star as well as its apsidal motion rate in a close binary system 
vary with the fundamental parameters such as the mass ratio and orbital distance. More specifi- 
cally, how good are the approximate methods described above for highly distorted stars? What is 
the relative role of tidal force to centrifugal force in increasing the central density concentration? 
Can the classical theory of apsidal motions based on perturbation analyses be accurately applied 
even to critical configurations? To address these questions, we in this work consider a circular-orbit 
binary in which two components are modeled as a polytrope and a point-mass, respectively, and 
calculate the equilibrium structure of the polytrope using an SCF method. The binary is assumed 
to be in spin-orbit synchronization, so that the tide is independent of time. As explained below, 
two salient features of our SCF method compared to those in the previous studies are that the 
former considers only the polytropic component in density computation, while explicitly including 
the tidal force from its companion as a fixed gravitational potential, and that the orbital velocity 
indispensable for calculating the apsidal motion rates is now obtained self-consistently from the 
equilibrium condition. 

The organization of this paper is as follows. In ^ we describe our implementation of the SCF 
method to find the equilibrium structure of a synchronously rotating binary component. In §3, we 
first present the numerical solutions for unperturbed polytropes and uniformly-rotating polytropes 
and compare them with the published results. We then turn to the solutions for both rotationally- 
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and tidally-distorted polytropes, and give detailed quantitative results on the changes of internal 
structures as well as on the parameters describing critical configurations. Using the self-consistent 
solutions of perturbed polytropes, we test the classical theory of apsidal motions in §4. Finally in 
§5, we summarize our results and discuss the astronomical implication of our work. 



2. Formulation 

2.1. Basic Equations 

In this paper we use an SCF method to study the equilibrium structure of a self-gravitating 
star under the external gravitational force from a companion in a close binary system. We assume 
that two stars revolve in circular orbits with constant angular speed Q and orbital separation A, 
so that the tidal force is stationary. We further assume that a star in question with mass Mi 
and volume radius i?i (hereafter star I), composed of inviscid polytropic gas without magnetic 
fields, rotates rigidly and synchronously with the orbital motion. For simplicity, we treat the 
companion star (hereafter star II) as a point mass with mass M2. We work in a rotating Cartesian 
frame at 0, whose origin lies at the center of mass of star I. Figure [T] schematically illustrates the 
binary configuration and the coordinate system in which the x-axis points toward star II located 
at (x, y, z) = {A, 0, 0) and the z-axis is perpendicular to the orbital plane. Star I is distorted 
rotationally and tidally, and has a polar radius Rp and an equatorial radius Re along the positive 
X-axis. We assume for computational convenience that all the substance of star I is contained 
within a sphere with radius Rl, the distance to the inner Lagrange point from the origin. 

The equation of hydrostatic equilibrium then takes the form 

VP 

+ V$ = 0, (2) 

P 

where P is the pressure, p is the density, and <I> is the effective potential defined by 

= $n + ^* + ^g- (3) 
In equation dSj), $q is the centrifugal potential 

^n = -^{[x-A{l-p)f + y% (4) 

with fi = Ml /{Ml + M2) being the mass fraction of star I, <I>^, is the gravitational potential from 

star II ^ , ^ 

^ ^ GM2 

i/(x - + y2 + ^2 ' 
and is the self-gravitational potential of star I that satisfies the Poisson equation 



(6) 
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The proper boundary conditions on are that and its first derivatives should behave well at 
r = (x^ + + z'^y^'^ — > and r — > oo, and should be continuous on the surface of a polytrope. 
These boundary conditions are automatically satisfied when we use the integral representation of 
the Poisson equation in ^2.31 below. 

We assume that the matter of star I obeys a polytropic relation 

P = (7) 

with the pressure constant K and the polytropic index n. Equation ([2]) is then integrated to yield 

$ + (n + = (8) 

where is the potential on the surface of star I with vanishing density. An equilibrium condition 
for star II is simply 



dx 



(9) 



x=(A,0,0) 



which must be checked a posteriori. 



Following the convention, we introduce the dimensionless coordinates ^ = rj, (") = (x, y, z) /a, 
where 



a 



{n + l)K 
AirG 



1/2 



We also define the generalized Lane-Emden function Q through 
where the central density pc is given by 

Ml 

The equilibrium condition ([8]) then becomes 

6 = ^',-^', 

where ^I' = ^/{[n + l]Kpc^^) and is the value of on the surface of star I. 
In dimensionless forms, equations ([I])-® and Q read 



(10) 



(11) 



(12) 



(13) 



(14) 
(15) 
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and 



9", 



(16) 



and (US]), nj^ 



where V denotes a dimensionless derivative with respect to ^. In equations ([T] 
G{Mi + M2)/A^ is the Keplerian angular velocity, u = A/a is the dimensionless orbital distance, 
and V is the dimensionless angular velocity defined by 



2ttGpc 




(17) 



which is a measure of the amplitude of rotational perturbations (e.g., IChandrasekhaiill933bl ). It is 
sometimes useful to measure the degree of rotational perturbations in terms of the mean density 
p = 3Mi/(47ri?f) as 

02 _ 2n^Rl 



X 



2ttGp 3GMi 



(18) 



(e.g., IStotherslll974l : IClaretlll999l l. 



Because of tidal and rotational distortions, the orbital motions of close binary stars are different 
from purely Keplerian ones. For binaries in spin-orbit synchronization, the angular velocity is 
naturally constrained by equation Q, or its dimensionless form 




(19) 



4={i.,0,0) 



However, it turned out in our SCF method that calculated from equation (|19|) was too small 

and varied too much to give converged results at the very early stage of iterations. Instead, we 
empirically found that the condition 




1 + 



1 



1/(1-^) 9e 



(20) 



corresponding to a vanishingly-small pressure gradient at the center of mass of star I, allows the 
solutions to converge rapidly, and is in fact very similar to equation (jl9p when an equilibrium 
is achieved. In practice, therefore, we first use equation (|20p to update until a temporary 
convergence is made, and then switch to equation (jl9p to obtain the desired equilibrium solutions. 

To check the accuracy of equilibrium configurations that we construct, we use the virial pa- 
rameter e defined by 



2T + 3jPdV + W 



W 



(21) 



where T = iO^Ma^V^ - / P^nd^ X is the total kinetic energy and W = J ($g -|- ^:^)pd^x is the 
total gravitational potential energy. 
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2.2. Structure Parameters 

A polytrope distorted by tidal and centrifugal forces in general has different radii in the equa- 
torial and polar directions. To measure its size, we use the volume radius 

/ \ 1/3 

where V denotes the volume enclosed by the isodensity surface with = 0. Note that in the absence 
of perturbing forces, is equal to the first zero of the Lane-Emden functions for unperturbed 



polytropes with n < 5 (e.g., IChandrasekhaiill939l ). 



Substituting equation (jl2p into equation (jlOp and eliminating a in the resulting equation by 
use of equation (j22p yield the mass-radius relation 

K = 7V„GMf -'^/"i^f-")/", (23) 

where Nn denotes the numerical coefficient 

{n-3)/n / \ (1-)/" 

= r , / e^d^e . (24) 



(n+ 1 

We note that equation (j23p holds also for unperturbed polytropes, but with different values of Nn 
To characterize the degree of the mass concentration toward the center, we use 

Pc 47r,ei 



p 3/e"d3e' 

The dimensionless moment of inertia of star I about the z-axis is defined by 

1 



(25) 



a^Pc 



p{x^ + y^)d'^x. (26) 



Finally, the total angular momentum of star I is given by 



2 



^1 = j P^ud-'x = Jorb + Jspin, (27) 

where Jorb = OMiA^(l — /i)^ and Jspin = ^ f p{x^ + y'^)d^x = a^pc^I are the orbital and spin 
angular momenta of star I, respectively. We will study in the next section the dependences of ^i, 
^n, Pc/Pi I, and Jspin/ Jorb upon the strength of centrifugal and tidal forces. 



2.3. Expansion of Density-Potential Pair 



For given values of n, p, and z^, equations (fT3P - (fT6]) constitute a closed set of equations for the 
equilibrium density 0"; the rotation parameter v is related to the dimensionless angular velocity 
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through equation ()17p and Q/i^x is constrained by equation (jl9p for binaries in spin-orbit 
synchronization. In an SCF method, the density and the potential are updated alternately and 
iteratively in such a way that at the i-th stage of iteration, equations (fn|) -p ^ are used to obtain 
the potential from Q,;, which is th en substituted in equation (jl3p to yield ©j+i for the next 
iteration stage (jOstriker &: Marklll968l ). 

We solve equation (jl6p using a mulitpole expansion technique in the region of space bounded 
by a sphere with radius Rl- For this, it is convenient to set up spherical coordinates {r,9,ip) 
centered at the center of mass of star I. We divide the sphere into N^i concentric shells and expand 
the density using the associated Legendre functions P^(cos0) as 



oo k 



(28) 



k=0 1=0 

where pki is the coefficients to be determined, and 

Rl 

' aNR 



(j - 1/2) 



(29) 



denotes the radius of the j-th shell for j = 1, 2, • • • , Nji. Note that equation (j28p does not contain 
sml<j) terms because of the even symmetry about the meridional plane. Using the orthogonal 
properties of the Legendre functions, one can show that 



Pkl[rj, 



2k + I 
47r 



TT 277 



i'^-^io)^^J J e''{rj,e,ip)pUcose)cosi^smededv, (so) 





where 5io is the Kronecker delta: Sio = 1 for / = 0; 5io = otherwise (see e.g.. iHernquist fc Ostriker 



19921 ^. 

For 0" expanded in equation (j28p . equation (jl6p yields the series solution 

oo k 



kl[r)r 



-k-l 



+ Rki{r)r^ pI.{cos 9) cos lip, 



(31) 



fc=0 1=0 



where 



Qklir) 



1 



2k + 1 



r'^^^Pki{T')dr', 



(32) 



and 



Rkl{r) 



1 



2k + I 



r'^-^PklirVr', 



(33) 
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(see e.g., iBinnev &: Tremaind l2008l ) . Note that equation (j3ip gives the gravitational potential at 
not only the interior but also exterior of a polytrope; a t r > Rr, can be obtain ed by taking 
Qki{r) = QkiifNu) and Rkiir) = RkiirNn) = 0. As noted bv lOstriker &: Mara (|l968l ). the integral 
representation via equations (p8]) - ([33]) of the Poisson equation is favored over the differential 
form (eq. [16]) since the former not only incorporates the boundary conditions but also guarantees 
numerical convergence. 



2.4. Iteration Procedure 

To find the self-consistent equilibrium solution for a given set of the parameters n, /i, and 
I/, we proceed the following steps. First, as a trial density distribution, we consider a uniform 
sphere with radius Rl/o- and density B = 1, consisting of = 1200 shells. The sphere is initially 
rotating at = 1- Second, we evaluate t he integrals in eq uations ([30|) . ([32|) and ([33]) using 



Gaussian and Newton-Cotes quadratures (e.g.. |Press et al. Ill988l ) to obtain the coefficients Pkii^j), 
Qikh and Riki successively. The corresponding self-gravitational potential is then obtained from 



equation ()3ip in which the summation over the /c-index is truncated after 10 terms; the associated 
error is negligibly small. Third, we calculate the rotation parameter v from equation ()17p and 
construct the effective potential ^. Fourth, we solve equations (fT3]) and (f20l) (or, equation [19] 
after a temporary convergence) to update G and fi/ilx- Fifth, we calculate the center of mass 
^cm = (/ / {J Q^d?i) of star I using the updated 0, and shift the coordinates so as to make 

coincide with the origin, while keeping star II at ^ = (i^, 0, 0) in the new coordinates. As a 
convergence criterion, we require the virial parameter e from equation (j2ip to be less than 10^^^; 
otherwise, we go back to the second step with the updated O and ^/VLk in the shifted coordinates, 
and repeat the iterations until the convergence is attained. For stars with critical rotation, we find 
that typically 50 iterations suffice to lead to converged results. 



3. Numerical Results 

While the main purpose of this work is to find the changes in the internal structure of binary 
stars due to rotational and tidal perturbations, we first apply our method to construct equilibrium 
models of polytropes in isolation or in uniform rotation. This will allow us to check our SCF method 
as well as to quantify the effects of rotation on the stellar structure. We then explore equilibrium 
polytropes subject to both centrifugal and tidal forces. 



3.1. Undisturbed Polytropes 



We first consider undisturbed polytropes, the solutions to which can be obtained by taking 
u = 0, corresponding to = = 0, in our SCF method. Table[T]lists the resulting dimensionless 
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parameters that w e find for some selecte d values of n. All the values are in excellent agreement 
with those given in IChandrasekhaii (|l939l ). confirming the performance of our technique. 



3.2. Uniformly Rotating Polytropes 



If one takes /U = 1 in equations p4j) and ([15]) and regards f as a free parameter, the problem 
is reduced to finding solutions for polytropes under uniform rotation. Figures [2] and [3] give as solid 
lines the dependences upon the rotation parameter v of the mass concentration pc/ p, dimensionless 
volume radius ^i, polar radius = Rp/a, and equatorial radius = Re/ct for rotating polytropes 
with n = 3 and 1.5, respectively. For comparison, w e plot the results availa ble from the published 
work, based on the first-o rder perturbation theory (jChandrasekhaiill933al). the second -order ap- 

and the 



proximation (jAnandlll968l ). the a dvanced first -order approximation (jLinnell 



1977 



1981 



Legendre-polynomial expansions (jjameslll964l ) as dotted lines, dot-dashed lines, open circles, and 
filled circles, respectively. 

As V increases from zero, the centrifugal force tends to distort the shape of a polytrope by 
increasing the equatorial radius at the expense of the polar radius. The polar compression is smaller 
than the equatorial expan sion, causing the polytrope to have a larger mean radius than the non- 
rotating counterpart (e.g.. 



Linnel] 



1981 



). A polytrope with faster rotation should definitely be more 
centrally condensed in order for self-gravity to offset the enhanced centrifugal and pressure gradient 
forces near the equatorial plane, as Figures [2]a and[3]a evidence. When v becomes sufficiently large, 
perturbed polytropes reach a critical state where they are on the verge of equatorial breakup. This 
critical rotation is formally expressed as d^/dr = at the equator, or equivalently, 96/5^1^=^^ = 
when V = Vc, indicating that the centrifugal force barely balances the self-gravity at the equator. No 
equilibrium configuration exists for v > v^- Table [2] gives the various parameters for the critically- 
distorted uniformly-rotating polytropes with differing n. The smaller n, the smaller Pc/p and the 
larger deviation of the gravitational potential of star I from that of a Newtonian point-mass. The 
values of Ac for polyt ropes with ind ices n = 1 .0—2.5 agree, within less than 0.2%, with the published 
numerical results of I JamesI (jl962l ) (see also iHurlev &: Robertslll964l ). Note that Ac = 0.3 6074 for 
n = 4.5 polytropes are virtually the same as the value from the Roche model with p = I (|Horedt 



20041 ) 



Although the first- and second-order perturbation methods are reasonably good for since 
the centrifugal force is very small near the poles, they systematically underestimate Pc/p, (,e, and 
^1 as V increases toward Vc- Being valid only for v <^ 1, they are also unable to yield correct values 
of Vc- On the other hand, f or the the wh ole range of v, our results for ^p and are within 1% 
of the numerical values from I Jamea (|1964| ). verifying again the accuracy of our SCF method. For 
practical purposes, we fit the dependences upon v/vc of the various parameters characterizing the 
shapes and internal structures of rotating polytropes using a simple function 



f{v) = ao + ai(^^^+a2(^^^ +a3ln(l 



04 



V 
LVr. 



(34) 
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with coefficients Oj. Table [3] lists the fitting coefficients for f = Ci? ^e, Cp, Pc/p, ^n, and / of rotating 
polytropes with n = 1.5 and n = 3; the associated fitting errors are less than 1%. Note that qq in 
each fit is identical to the corresponding value of non-rotating polytropes. 

It is interesting to compare the shapes of rotationally-distorted polytropes with those under 
the Roche approximation that assumes that the stellar gravity is dominated by the central mass 
concentration and that the polar radius is unchanged due to rotation. The main predictions of the 
Roche models are that the ratio of the equatorial velocity Ve to the critical value Vc behaves with 
the aspect ratio TZ = Ce/S,p as 



Vr 



3(7^- 1) 
n 



(35) 



and that TZr = 1-5 iii the case of critical rotation (e.g., ICollind 119631 : iBaumgarte k. Shapiro! 11993 : 



Townsend et al.ll2004l ). Figure |l]a plots as a solid line the change in TZc of rotating polytropes with 
varying n, while Figure gives Ve/Vc as a function of TZ for various polytropic and the Roche 
models. In terms of TZc, the Roche models are a good approximation for polytropes with very high 
mass concentration toward their centers (relative errors less than 3% for n ^ 2.0). On the other 
hand, equation (i35]l is quite accurate for polytropes with n = 1.5 except near the critical config- 
uration. It underestimates Ve/Vc for polytropes with n ^ 3.0 by about 10%, while overestimating 
for n = 0.5 polytropes by about 15%. This suggests that one should be cautious when estimating 
the equatorial velocity from an observed aspect ratio of a rapidly rotating star. 



3.3. Rotationally and Tidally Distorted Polytropes 

When a binary star with a given polytropic index n is subject to both tidal and centrifugal 
distortions, its internal structure is determined solely by the two parameters in our formulation: 
its relative mass p and the dimensionless orbital separation v; for systems in synchronous rotation, 
the rotation parameter v is given automatically via equations ()17|) and ()19p . Figure [5] plots the 
dependences of v upon v for the equilibrium polytropes under both tidal and centrifugal forces. Each 
solid line corresponds to a sequence of equilibrium models with a fixed relative mass p displayed 
in the panels. Note that v monotonically decreases with increasing orbital separation along each 
sequence, becoming vanishingly small at — > oo. On the other hand, a polytrope becomes distorted 
more and more as v decreases along a given sequence. When v becomes sufficiently small, it 
eventually reaches a critical state beyond which self-gravity no longer balances the combined tidal 
and centrifugal forces at the equator. The filled circles marked at the uppermost tips of the 
sequences correspond to the critical points (ucVc)- In what follows, we use the term "critical 
component" to refer to star I in critical configuration. The dotted line connecting the critical 
points draws the locus of the critical components in the v-v plane, above which no equilibrium 
polytrope exists. 
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3.3.1. Shape, Size, and Internal Structure 

Figure [6]a-c illustrates how the shapes on the meridional plane of n = 3 equilibrium polytropes 
change as the amplitude of tidal and rotational perturbations varies. All the models have ^ = 0.5. 
The corresponding Roche lobe around star I is plotted as a dotted line. In each panel, the inner 
solid contours draw equidensity surfaces with p/p = 3~'"(pc/p)iso! where m = 0, 1, • • • ,5 from 
inside to outside, while the outermost contour represents the boundary of a perturbed polytrope. 
Here and in what follows, the subscript "iso" indicates the quantities of the corresponding isolated 
polytrope. When v is sufficiently large, the tidal and rotational perturbations are so weak that a 
polytrope is in an almost spherical configuration, with the central density only slightly larger than 
the unperturbed value (Fig. [6]a). As v decreases, the enhanced perturbations make the polytrope 
larger and more centrally concentrated (Fig. [6]&). It is remarkable that when u = Vc {= 18.81 for 
n = 3), the outer envelope of the critical component fills the Roche lobe almost completely, while 
the dense inner-part still remains nearly spherical (Fig. [6]c). 

Figure [Hd plots the outer boundaries of the critical components with differing n for /i = 0.5. 
When critical component has a boundary virtually identical to its Roche lobe. With 

relatively low central concentrations, polytropes with n ^ 1 become flattened toward the equatorial 
plane, although the deviation from the Roche lobe is still slight. Figure [7] gives a quantitative 
comparison between the volume radius of a critical component and the effective radius Rl 



the Roche lobe fitted bv lEggletonI (11983) 



Rl ^ 0.49g^/3 

A 0.6(?2/3 + ln(l + (;V3)' ^ ^ 

where q = M1/M2 = /u/(l — p). Indeed, Eggleton's formula is an excellent approximation to Ri, 
with a relative error less than 2%, for polytropes with n ^ 2. For small values of n, Rl overestimates 
Rl because of the equatorial flattening. For n = 3/2 polytropes, we find 

Rl 0.5126g°-^388 _ 

A ~ 0.6710(?0-7349 + ln(l + gO.3983) ' ^o"^ n- 61 

gives a good fit (within ~ 1 — 2% for 10~^ < Q < 10^) to the volume radius. 

Figures El and [9] plot the variations of Pc/p, ^1, Nn, and / of rotationally- and tidally-distorted 
polytropes with n = 1.5 and 3.0 against the relative radius Ra = Ci/^ = Ri/A as solid lines. Again, 
the dotted line connecting the top ends of the constant-/^ sequences in each panel corresponds to the 
critical components. As expected, the central density concentration, volume radius, and moment 
of inertia increase as the perturbation amplitude increases. The relative increment of Pc/p is 
greater for polytropes with larger n. Note that the largest changes in these quantities occur when 
p = I, corresponding to purely rotating polytropes; in synchronized binaries with 0.1 ^ p ^ 0.9, the 
variations of Pc/p, Cii are less than 10% even in critical configurations. The dependences 

of these quantities on Ra for a fixed value of p can be approximated by 

g{RA) = 60 + hiRA + b2R\ + hR\ exp (64^?^), (38) 
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with g = Pc/p, Ci) or /. The fitting coefficients hi are given in Tables HI O [6l andO respectively; 
the fits are accurate within 1%. Note that 5o's in each Table are identical to the values for the 
unperturbed polytropes. 

Since a depends on Pq which in turn varies with v and /x, the dimensionless volume radius 
shown in Figure [8] does not give direct information on the change of a stellar size due to centrifugal 
and tidal forces. Assuming that the perturba tions do not modify the p ressure constant K (equiva- 



lently, specific entropy or degeneracy; see e.g.. iBurrows &: Liebertlll993l ) inside a star, equation ([23 



gives iii/i?i,iso = {Nn/Nn^iso)"'^^^~^^ as the size ratio of distorted to unperturbed polytopes with 
the same mass. Figure [TO] plots Ri/Ri^iso as solid lines for polytrope with n = 1.5 and 3.5, showing 
that stars become bigger as the amplitude of perturbations increases^, For n = 3.5 polytropes in 
critical configuration, the enhancement factor in the stellar size is ~ 12% for p = 1 and ~ 2 — 3% 
for p ~ 0.1 — 0.9, while it increases to ~ 23% for p = 1 and ~ 4 — 8% for p ~ 0.1 — 0.9 for n = 1.5 
polytropes. Taken together with the results for pc/p discussed above, this implies that polytropes 
with smaller n adjust themselves to the perturbations by increasing their size more than being 
centrally condensed, while those with larger n do so by increasing the density concentration more. 

In the dynamical study of close binary systems such as mass transfer, the spin angular momen- 
tum is often neglected compared to the orbital angular momentum. To check if this is a reasonable 
assumption, we plot in FiguredUthe dependences of Jspin/^orb on and p as solid lines. Overall, 
J spin/ Jorh IS Well approximated by 



«^spin 
Joih 




(39) 



with the coefficient cq = 0.20460 and 0.07536 for n = 1.5 and 3 cases, respectivelyl^ The fit is 
almost exact for small and the maximum error occurring at the critical values of ^i/z^ is less 

than 5%. Apparently, the contribution of the spin angular momentum to the total is negligible for 
wide binaries where ^i/j^ <^ 1. Notice, however, that when a binary star with sufficiently large p 
is near the critical configuration, Jspin even exceed Jorb and can thus be dynamically important. 

For the critical components, Pc/p, Nn, Jspin/Jovh-, and / can be approximated within 1% of 
accuracy by 

0.27555°-^^^^ 

Pc/p = 5.9907 H „ ; — ; TTTTET, for n = 1.5, (40a) 

^'^ 0.2430g0-8258 +in(l + gi-i4i5)' V ; 

5 3141g°-^^^'5 

pj-p = 54.1825 + 3,33^^^^^^^,,,,,^ , forn = 3, (40b) 



^ As is well know n, it is not viable to constrain the radius from the mass for polytropes with n 
Chandrasekhar 19391 ). 



= 3 (e.j 



^In the limit of ^i/;^ — > 0, it can be shown that co = 3J(pc/p)(47r5f ) ^ for unperturbed polytropes. 
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and 



0191o°-^^^^ 

^" = 0.3557,0.9130 + 1,(1 + ,i.ia35) ^ ^orn = 1.5, (41a) 

= °-^^^^'- 0.35430gi-00296+ 1,(1 + ^1.08931) ^ f-- = 3, (41b) 



71-858 + 0.043g0-645^ 

-^orb \o.068(g - 0.5)2-i°6 + 0.100g°•^ if g > O.5., 

^spin _ f0.026gi-853 + 0.015g°-^5\ if g < 0.5, 

^rb \o.026(g - 0.5)2 °6i + 0.037g°•^ if g > 0.5, 



Jspin J0.076gi-858 + 0.043gO-645^ if g < 0.5, 

— for n = 1.5, (42a) 



for n = 3, (42b) 



15.5476g0-9334 

^ = 0.3795,0.8902 + 1,(1 + ,1.1985)' fom = 1.5, (43a) 

4 4373a°-^8^^ 

I = 90.910 + , , 7T^, for n = 3, (43b) 

0.5200,0.9560 + 1,(1 + ,1.1593 ' ' V ; 



for the range of 0.05 < , < 40. 



3.3.2. Rotational vs. Tidal Distortion 



Figure [T2] plots as solid curves Vc and Uc for the critical components with n = 1.5 and 3.0 as 
functions of fx. These are fitted to 

^3 - 4.08028 + 4.10(1 - ;u) - 0.54(1 - /x)2 - 6.58(1 - /i)°-^^ for n = 1.5, (44a) 
: 4.35828 + 4.55(1 -/Li)- 0.79(1 -//)2- 7.15(1 - //)°-^3 forn = 3, (44b) 



Vc 

Vn 



10^ 



and 



T.l=od^{l- ifo.i</x<i, 

-1/3 



(45) 

eo + ein ^'-^ if ^ < 0.1, 

wi th the coefficients dj and Ci given in Table El Filled circles in Figure [T2] a represent the results 
of ISingh &: Singhl (|l983l ) un der the double a pprox imation, which are within ^ 10% of our fully 
self-consistent values. While ISingh &: Singhl (jl983l ) found that Vc monotonically increases with 



( ^ 0.09), our results show that Vc attains a minimum at /ii. 



0.13-0.14forn > 1.5. 



For /i > /Xmin, the rotational eff ect is inore im portant than the tidal force in maintaining 
equilibria of critical components (e.g., lRussel]||l928l ) rl As /u increases, star I would feel a weaker 



iRusselll (|l928l ) considered the ratio of the tidal to centrifugal forces at the undistorted stellar surface, and found 
that tidal and rotational effects are the same when p = 1/6, which is not much different from /imin that we have 
found. 
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tidal force from its companion if the orbital distance were the same. Therefore, star I must rotate 
faster to be critical, increasing Vc- This in turn corresponds to decreasing i/c due to the condition of 
spin-orbit synchronization. In the limit of — > 1, converges to the value found in ^3.2l for critical 
rotation in the absence of the tidal force, while i^c is tending to suggesting that a fictitious 
secondary with vanishing mass is located at the equator of the primary in critical condition. 

When fi < /imin; on the other hand, the tidal force dominates the centrifugal force in the force 
balance. Note that v is proportional to the ratio of the tidal force to the self-gravity at r ~ so 
that Vc oc {S^i/ucY / ^. Eggleton's formula (eg. [36]) shows that ^i/z^c decreases slightly more slowly 
than /i^/^ as ^ decreases. This causes Vc to increases slowly as /i decreases from /imin- In the limit of 
^ ^ 0, z^c oc and Vc converges to 1.1530 x 10~^ and 1.2496 x 10"'^ for the cases with n = 1.5 

and 3.0, respectivelylll Figure [13] plots the variation of /Umin with the polytropic index, which shows 
that fXrain = 0.136 as n approaches 50 For lower values of n, the enhanced equatorial radii due 
to flattening (e.g., Fig. [6|(i) make the centrifugal force more important at the surface, tending to 
decrease ^min- 



3.3.3. Structure Constant k2 

Once the three-dimensional density distributions of disturbed polytropes are found, it is 
straightforward to calculate the apsidal motion constant k2 defined by 



3-r?2(gi) 
4 + 2j?2(6) 



(46) 



where ??2(Ci) is the solution at the volume radius of the Radau equation th at describes the distortion 
of equidensity surfaces inside a perturbed polytrope (see, e.g.. lKopallll978l ). For isolated polytropes, 
we have confirmed that equation (|16]) yields log A;2,iso = —0.8438, —1.8403, and —2.3081 for n = 1.5, 
3.0, and 3.5 cases, r espectively, which are in good agreement with the published results (e.g., 
Brooker fc 01l3ll95d l. 



Figure [131 plots as solid lines the changes of A log k2 = log(A;2/A;2,iso) with respect to the rotation 
parameter A and the central density concentration Alog pc/p = log[{pc/ p) / {pc/ p)iso] for distorted 
polytropes with n = 1.5, 3.0, and 3.5. Clearly, /c2 decreases as the perturbation amplitude (or, 
the degree of central concentration) increases. Using realistic ZAMS models with modern opacities 
for unperturbed stars, but by relying on a quasi-spherical or double-approximation method in 



*For the Roche model, Horedtl 1 2004 ) showed that the mean critical rotation parameter is Ac ~ 2*/3^ for a 



vanishingly small (see also Kopall 



19891). This corresponds to Vc = Xc/{pc/p) ~ 1.002 x 10"^ and 1.041 x 10" 



since pc/p ~ 6.5721 and 63.2685 when fi —> for n — 1.5 and 3.0, respectively. These are in reasonable agreement 
with our results for the critical components. 

^Using the Roche model, Kopall ( 19781 ) found the minimum value of Ac occurs in between p, = 0.130 and 0.167 (or 
g = 0.15 and 2), which is consistent with our results. 
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evaluating the sole effect of stellar rotation, IStothersI (|1974| ) and IClaretl (|l999l ) found A log k2 = 
— 0.7A and Alog A;2 = — 0.87A, respectively. These are plotted as dotted and dashed lines in the left 
panels of Figure [U approximately e nveloping our re sults for n = 3.5. While it is not straightforward 
to compare our results with those of lStothersI (jl974l ) and lClaretl (j 19991 ) since the unperturbed stellar 
models are different, this suggests that rotationally-distorted stellar models considered by these 
authors likely have internal structures similar to those of n = 3.5 polytropes. The right panels 
of Figure [14] show that for polytropes with fixed n, AlogA;2 varies approximately linearly with 
Alogpc/p, nearly independent of /i, with an average slope of —1.12, —1.33, and —1.39 for n = 1.5, 
3.0, and 3.5 cases, respective ly. This confir ms the notion that k2 is inversely proportional to the 
central density concentration (|Stotherslll974l ). Figure [HI also shows that |Alog A;2| < 0.1 for fj, ^ 0.9, 
while it can be as large as 0.3 for strongly centrally-condensed polytropes in critical rotation with 
fj, = 1. This suggests that the decrease in /c2 due to both rotational and tidal distortions is much 
smaller in binary stars with synchronous rotation than in purely rotating stars. 



4. Apsidal Motion 



Stars in close binaries have a shape that deviates from a spherical configuration, rendering 
the orbits of their companions non-Kepler ian. As mentioned in the Introduction, it is well known 
that tidal and rotational distortions cause the apsidal line to advance, the rate of which depends 
on the degr ee of ste l lar co ncent ration . Man y observational studies have employe d the analytic 
formulae of ICowlingI (|l938l ) and ISternI (|l939l ) , which are an eccentric extension of iRusselll (|l928l ) 
formula applicable to near-circular orbits. By retaining only the linear-order terms, however, all of 
these formulae are valid in a strict sense when the tidal and rotational distortions are very small. 
Even in the cases when the effects of dynamical tides and relativity are unimportant, however, the 
distortion of a stellar shape can be substantial e specially for cr itical components. In this section, 
we check the validity of the classical formula of iRusselll ()1928l ) using self-consistent solutions for 
stars in circular-orbit binaries for which tides are stationary. 

To obtain the apsidal motion rates from our distorted polytropes, we calculate the self- 
gravitational potential $g numerically along the axis joining the two binary components after 
the self-consistent solutions are obtained. We then apply the orbit theory in which the rate of 
change of the periastron is given by 

n-K, (47) 



where 



K = 



Ml A3 dx 



(48) 



=(A,0,0) 



is the epicycle frequency at the location of star II (e.g., iMurrav Dermottlll999l ). The resulting 
(numerical) apsidal motion rates, tn^, versus the relative radius Ri/A = ^i/v are shown as solid 
lines in Figure [15] for = 0.1, 1.0 and n = 1.5, 3.0, 3.5 cases. Also plotted as dotted lines are 
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the (theoretical) rates, wt, from the classical formula (eq. [T]) with k2 calculated in N3.3.3i When 
the relative radius is small (i.e., 5,i/v < 0.15 for /x = 0.1 and ^i/j^ < 0.5 for ^ = 1) and thus the 
rotational and tidal perturbations are weak, the agreement between zu„ and wt are remarkably 
good. As the amplitude of perturbations increases further, however, the classical formula begins to 
underestimate the real apsidal motion rates. This is because the the distortion of star I becomes 
no longer in the linear regime as ^i/i^ approaches the critical values. In addition, we found that 
the gravitational potential due to a highly distorted component has an non-negligible contribution 
fro m terms prop ortional to while the classical theory took only or higher-order terms (see, 



e.g. 



Sternlll939l ). making rij^ smaller than w„ 



Figure [T6l plots the relative difference {wn — wt)/wt of the numerical and theoretical rates of 
the apsidal motions for critical configurations as a function of jj,. For critically distorted polytropes 
with n = 1.5, the real apsidal motion rates are larger by ^ 50% than the values predicted by 
the classical theory. As n increases, the fractional mass occupying outer layers that exhibit large 
distortions decreases, yielding {vjn — wtj/wt ~ 5 — 15% for stars with n = 3.5. These results 
suggest that the classical formula still works quite well for stars with large density concentration at 
centers, while it underestimates the apsidal motion rates significantly for stars with relatively low 
Pc/p, such as in low-mass convective stars. 

Finally, we define the "effective" internal structure constant A;2,cfT through 

(49) 




We note that /i;2,efr introduced in equation (p9|) is simply to match the real apsidal motion rates by 
means of the classical theory with the modified values of A:2, and thus it should not be interpreted 
as conveying physical information on the internal structure of perturbed stars. Figure [T7] shows 
the dependence of Alogfc2,efr = log(A;2,efr/^2,iso) on A for n = 1.5 and 3.5 polytropes. Dotted lines 
correspond to the cases with critically distorted polytropes. Note that A log A;2^efr even increases 
with increasing A when n is small, although this does not mean stars become less centrally-condensed 
as the amplitude of perturbations increases. For practical uses, we fit A;2,efr using 

A log = ai A + aaA^ + agA^, (50) 

with the coefficients ai listed in Table [9] for n = 1.5, 3.0, and 3.5 cases. The relative errors of the 
fits are less than 2%. 



5. Summary and Discussion 

Binary stars are a traditional source of information on stellar mass, radius, luminosity, etc., the 
inter-relationships of which are of crucial importance to stellar astrophysics. Many observational 
and theoretical arguments suggest that centrifugal and tidal forces may significantly affect the 
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physica l properties of stars in close binary systems. Since the pioneering work of IChandrasekhar 
(|l933al ffl ) , there have been numerous theoretical studies that tried to find equilibrium solutions of 



polytropic stars under rotational and tidal perturbations. However, most of them take some type 
of approximations, and thus are valid only when the perturbations are weak. In order to construct 
equilibrium structures of such polytropes without relying on any assumption on the perturbation 
strength and the degree of stellar distortion, we use an SCF method in which the angular orbital 
velocity $7 is self-consistently constrained rather than being treated as a free parameter. 

To ignore the complicated effects of time-dependent tides and to make contact with the previous 
theoretical work, we consider a circular-orbit, synchronized binary consisting of a polytrope and a 
point-mass object. Based on our SCF scheme described in ^ we solve the dimensionless steady- 
state momentum equation (jl3p and Poisson equation ()16p in the rotating frame at O, alternately 
and iteratively, until the converged solutions are obtained. We first apply our method to polytropes 
in isolation and subsequently in uniform rotation, and calculate the density concentration pd p, the 
volumn radius ^i, the equ atorial radius £p^ , and the polar radius these are in good a greem ent 



with the accurate results of lChandrasekhan (119391 ) for unperturbed polytropes and lJamed ()l964l ) for 



rotating polytropes, confirming the performance of our SCF method. When the rotation parameter 
V (eq. [17]) is sufficiently large, a rotating polytrope reaches a critical state, beyond which no 
equilibrium is possible. Equation and Table [3] provide fitting formulae of ^i, ^p-, Pc/p, 
Nn, and the moment of inertia / for polytropes with n = 1.5 or 3. By comparing the axis ratio 
TZ = ie/ip and the velocity ratio Ve/Vc with the predictions from the conventional Roche model, 
we find that in terms of the critical value T^o the Roche model is a reasonable approximation 
for rapidly rotating stars with high central density concentration (i.e., when n ^ 2.0), although it 
underestimates Ve/Vc by about 10% when n ^ 3.0. 

We then calculate equilibrium solutions of polytropes subject to both centrifugal and tidal 
forces. The rotational and tidal perturbations in general make disturbed stars larger in size and 
more condensed toward their centers. For n = 1.5 polytropes in critical configurations, for example, 
the distorted stars are larger by ~ 23% and ~ 4 — 8% for = 1 and p ~ 0.1 — 0.9, respectively, than 
the unperturbed counterparts with the same specific entropy. Here, p is the mass of a perturbed 
star relative to the total mass of a binary system. For more centrally-concentrated polytropes with 
n = 3.5, the enlargement factor is reduced to ^ 12% and ~ 2 — 3% for p = 1 and p ~ 0.1 — 0.9, 
respectively. In determining the structure of a critical component, the rotational effect is more 
important than the tidal effect provided p > pmin ~ 0.13 — 0.14. As in purely rotationally-distorted 
cases, the shape and size of both rotationally- and tidally-perturbed polytropes are well described 
by the Roche models as long as n ^2. Equation ([38]) and Tables HI O [H and [7] give functional fits 
of Pc/p, Nn, and / on the mass radio p and the relative radius Ri/A = ^xjv, while equation 
()39p is for the fit of Jspin/'^rb) the ratio of spin to orbital angular momenta of star I. Equations 
(|40l) - (l45]) and Table ([8]) give the dependences on p of pdp-, Nn, Jspm/ Jorh, I, Vc and fc for critical 
configurations. 



Using self-consistent equilibrium solutions of rotationally- and tidally-distorted polytropes. 
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we calculate the internal structure constant k2 by integrating the Radau equation as well as the 
apsidal motion rate Wn based on the orbit theory. We find that A log ^2 depends almost linearly 
upon A log /?c//0, with an average slope of —1.12, —1.33, and —1.39, nearly independent of for 
n = 1.5, 3.0, and 3.5 polytrope s , resp ectiv ely. The exp ression AlogA;2 ~ —(0.7 — 0.87)A found 
for rotating stars from IStotherd (|1974| ) and IClaretl (|l999l ) appears to be valid only for stars with 
very high density concentration, corresponding to n ~ 3.5 . The co i npari son of Wn calculated from 
our SCF method with wt from the classical formula of iRusselll (jl928l ) based on the first-order 
perturbation analyses shows that the latter is valid only when the amplitude of perturbations is 
small. For the cases of critical configurations where tidal and centrifugal forces are rather strong, 
the classical theory underestimates the real apsidal motions rates by as much as ~ 50% for n = 1.5 
polytropes and ~ 5 — 15% for n = 3.5, indicating that the nonlinear effects are by no means 
negligible in augmenting the apsidal motion rates. To this end, we define the effective internal 
structure constant k2^cs that gives the real apsidal motion rates by means of the classical formula. 
Equation (j50p and Table [9] present the algebraic fits to A:2,eff- 

Observations of eclipsing binary stars indicate that B-type binary stars are about 20% smaller 
in size than the single stars of the same spect ral type , whil e A- an d F-type binary stars are larger 



than the corre sponding iso lated stars (e.g., iMalkov 



2003 



20071 ). For M-type stars with mass 



0.4 — 0.8 M0, iRibasI (j2006l ) reported that observed radii of eclipsing binary stars are larger by 



5 — 10% than those of the theoretical, isolated models. On the other hand, our self-consistent 
calculations show that centrifugal and tidal perturbations make binary stars bigger in size. These 
observational and theoretical results are qualitatively reconciled with each other if one considers 
the f act that early-type (late-type) binary stars rotate slower (faster) than the isolated counterparts 
(e.g., iMcNallyi Il965l : iMalkovl Il993l . 120031 ). Slow rotations of early- type stars in close binaries are 
apparently a consequence of the spin-orbit synchronization occurring during their main-seque nce 



phase, induced by tidal torque d ue to either radiatiy e damping of dynamical tides (IZahnlll977l ) or 
large-scale mechanical currents (jTassoullll987l . Il988l ). In addition, early- type binary stars appear 
to have a low initial rotational velocity amountin g to ^ 130 km s~^, compared to ~ 200 km s~^ 



np 

for single stars, when they enter the ZAMS stage (jKhaliullin &: Khaliullinall2007l ). 



By analyzing the nea r-IR, long-baseline interferometri c data of the Be star Achernar which is 
known as a rapid rotator, iDomiciano de Souza et al.l (120031 ) reported that the star has a projected 
aspect ratio oi TZ = 1.56 ± 0.05. Since the projection always tends to decrease the measured 
aspect ratio, this implied that the real oblateness must be larger than the critical value under the 
Roche approximation. If Achernar can be modeled by a single, rigidly-rotating polytrope, Figure 
Ula suggests that it must have n < 2.0, which is unlikely given tha t Be stars may have a larger 
central density concentration than that of a polytrope with n = 3. I Jackson et al.l (|2004l ) showed 



that differentially-rotating stellar models with a more realistic equation of state allow equilibrium 
stars with TZ l a rger t han 1.5, compatible with the interferometric observations. On the other hand. 



Carciofi et al.l (j2008l ) found that the data are also consistent with a critically-rotating central star 



surrounded by a small equatorial disk produced by the ejected material from the central star. 
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If one component (star I) in a semidetached binary overfills its Roche lobe, it may rapidly 
lose mass to its companion in a few orbital periods. Such dynamical mass transfer is likely to 
create a common envelope bi nary, important for the evolution o f cataclysmic binaries and double 
degenerate white dwarfs (e.g., iPacz vhskil 1 1 9 76l : iFrank et al.ll2002l ). In finding the condition for the 
onset of the dynamical mass transfer, most previous work considered the conservation of the orbital 
angular momentum, neglecting the contribution of the spin angular momentum. Using t he Eggleton 
formula (eq. [36]) for the Roche lobe radius, for instance, iHjellming &: WebbinkI (jl987l ) found that 
the rapid conservative mass transfer from an n = 1.5 polytro pe to its cq r npact companion occurs 
if q = M1/M2 > Qc = 0.634, while t he use of the less accur ate Paczvhski ( |l97l|) formula, RjJA 



0.4622(g/[l + g])V3^ yields Qc = 2/3 (IRaPDaport et al.lll982l : lO'Souza et al.ll2006l : iMotl et al.ll2007l 1. 
We revisit this issue using the conservation of the total angular momentum that includes Jgpin given 



by equation (j39p . and find that the critical value of q for the dynamical mass transfer is increased 
to Qc = 0.693 and 0.736 under the Eggleton and Paczyhski formulae, respectively, for n = 1.5 
polytropes. Note that the value of Qc based on the Eggleton formula including Jgpin is similar to 
that under the Paczyhski formula without the effect of Jspin- Although the changes in Qc due to the 
inclusion of the spin angular momentum are not great since Jspin/ Jorb ~ 0.1 for g < 1 (or ^ < 0.5), 
the effect of Jspin can be significant in the processes of mass transfer if the donor-to-accretor mass 
ratio is large. 

Measuring apsidal motion rates is probably a unique means to probe the internal density con- 
centration of bin ary stars. For systems where the effects of general relativity (jGimenezj Il985l ) or a 
third body fe.g.. iBozkurt &: Degirmencill2007l ^ are unimportant, most previous studies adopted the 
classical theory of ISternI (jl939l ) with k2 calculated from isolated polytropes arid found that binar y 



stars are muc h more centrally condense d tha n predicted (e.g., 



More recently, IClaret &: Gimenez I (119931 ) and IClaret &: Willema (120031 ) showed the theoretical re 



Schwarzschildl Il958l : iKopall Il978l ) 



suits are in good agreement with observations if the changes of k2 due to rotation as well as the 
effects of the initial chemical composition and convective core overshooting are taken into account. 
However, the reduction of k2 considered in these papers amounts to AlogA;2 ~ —(0.2 — 0.3), which 
may be possible only for purely rotating cases with /i = 1. Even for highly centrally-concentrated 
polytropes with n = 3.5, Figure [T^ shows that |AlogA;2| < 0.1 for synchronous binaries in the mass 
range of 0.1 ^ ^ 0.9. In addition, our SCF analyses show that the classical formula derived under 
the linear approximation is likely to underestimate the real apsidal motion rates for stars with low 
central concentration, especially in critical configurations. All of these suggest that special care 
should be taken in interpreting observed data for apsidal motions of close binary systems where 
stellar rotation is slow. 
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Table 1. Parameters of undisturbed polytropes 



n 


^1 


Pc/P 




I 


0.5 


2.7527 


1.8352 


2.5235 


117.575 


1.0 


3.1416 


3.2899 


0.6366 


101.844 


1.5 


3.6538 


5.9907 


0.4242 


93.156 


2.0 


4.3529 


11.4025 


0.3648 


88.895 


2.5 


5.3553 


23.4065 


0.3515 


88.190 


3.0 


6.8969 


54.1825 


0.3639 


90.910 


3.5 


9.5358 


152.8837 


0.4010 


98.398 


4.0 


14.972 


622.4079 


0.4772 


114.275 


4.5 


31.837 


6189.4719 


0.6580 


152.617 



Note. — n: polytropic index; ^i: dimcnsionlcss 
volume radius; pc/p- central density concentra- 
tion; Nn- coefficient of the mass-radius relation; 
/: dimensionless moment of inertia 



Tabic 2. Parameters of uniformly-rotating polytropes in critical configuration 



n 






A,. 


A'„ 


/ 


l>r/P 


si 






0.5 


1.4409x10" 


-1 


0.29386 


1.42922 


587.284 


2.0394 


3.4754 


2.1779 


4.8489 


1.0 


8.3880x10" 


-2 


0.33967 


0.45025 


236.366 


4.0494 


3.7367 


2.6922 


4.8278 


1.5 


4.3622x10" 


-2 


0.35168 


0.34330 


153.088 


8.0619 


4.2676 


3.2962 


5.3526 


2.0 


2.1576x10" 


-2 


0.35694 


0.32178 


120.783 


16.5434 


5.0863 


4.0556 


6.2961 


2.5 


9.9305x10" 


-3 


0.35922 


0.32965 


106.637 


36.1737 


6.3014 


5.1000 


7.7549 


3.0 


4.0803x10" 


-3 


0.36019 


0.35728 


102.111 


88.2755 


8.1907 


6.6738 


10.0548 


3.5 


1.3861x10" 


-3 


0.36057 


0.40779 


105.249 


260.1278 


11.4334 


9.3419 


14.0214 


4.0 


3,2913x10" 


-4 


0.36071 


0.49864 


118.369 


1095.9536 


18.1078 


14.8100 


22.1991 


4.5 


3.2295x10" 


-6 


0.36074 


0.70213 


154.862 


11170.2095 


38.7737 


31.7202 


47.5299 



Note. — n: polytropic index; Vc &c \c critical values of the rotation parameters; Nn- coefficient 
of the mass-radius relation; I: dimensionless moment of inertia; Pc/p- central density concentration; 
^1, ^p, & dimensionless volume, polar, and equatorial radii, respectively 
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Table 3 


. Fittini 


z, coefficients for tlie 


various 


parameters of rotating 


z, polytropes 




/ 






n = 1.5 










n = 3 






ao 


ai 




0-3 


a4 


ao 


ai 


"2 




a4 




3.65375 


0.03588 


0.00452 


-0.11724 


0.21911 


6.89685 


0.30538 


0.18974 


-0.24088 


0.96237 




3.65375 


0.06834 


0.01192 


-0.25324 


0.22619 


6.89685 


0.40291 


0.47814 


-0.45359 


0.99175 


?P 


3.65375 - 


-0.07555 


-0.00142 








6.89685 


-0.22243 


-0.00057 








Pc/P 


5.99070 


0.01990 


0.01092 


-0.52212 


0.22119 


54.18248 


3.83480 


4.55374 


-7.62154 


0.96447 


N„ 


0.42422 - 


-0.00904 


-0.00057 


0.00996 


0.21821 


0.36394 


-0.00599 


-0.00066 








I 


93.15635 


2.05830 


0.25902 


-30.87279 


0.17755 


90.90975 


8.53385 


2.62075 









Note. — ai's are defined through equation ((34). 



Table 4. Fitting coefficients for pc/p of distorted polytropes 



p 




n = 1.5, bo 


= 5.99070 






n = 3, 6o = 


54.18248 




bi 


b2 




b4 


bi 


b2 


bs 


b4 


0.1 


-0.17092 


3.49535 


0.06614 


28.20734 


-2.46052 


59.12506 


9.44135 


17.55872 


0.2 


-0.12860 


2.13006 


0.03425 


21.79792 


-1.85425 


36.09363 


5.38216 


13.12862 


0.3 


-0.10830 


1.58568 


0.02566 


18.13797 


-1.55648 


26.79729 


4.29253 


10.56543 


0.4 


-0.09552 


1.28070 


0.02340 


15.36396 


-1.36230 


21.48995 


4.03468 


8.63850 


0.5 


-0.08637 


1.08067 


0.02478 


12.96732 


-1.21555 


17.90157 


4.23624 


7.02247 


0.6 


-0.07921 


0.93547 


0.03011 


10.74527 


-1.09285 


15.18838 


4.84080 


5.60038 


0.7 


-0.07320 


0.82138 


0.04173 


8.62004 


-0.98314 


12.95893 


5.89126 


4.33881 


0.8 


-0.06148 


0.68295 


0.08494 


6.11931 


-0.88629 


11.05241 


7.38448 


3.25984 


0.9 


-0.06648 


0.65962 


0.09659 


4.98583 


-0.83960 


9.68781 


8.75499 


2.46508 



Note. — bi's are defined through equation II38I I. 
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Table 5. Fitting coefficients for the volume radius of distorted polytropes 



M 




n = 1.5, bo 


= 3.65375 






n = 3, bo 


= 6.89685 




bi 






b4 


bi 


b2 


ba 


b4 


0.1 


-0.05174 


1.22025 


0.13999 


19.60622 


-0.10449 


2.60661 


0.64437 


15.45073 


0.2 


-0.03954 


0.75031 


0.07579 


14.87539 


-0.07951 


1.59668 


0.36468 


11.50900 


0.3 


-0.03346 


0.56039 


0.05810 


12.14228 


-0.06669 


1.18356 


0.29108 


9.18447 


0.4 


-0.02956 


0.45272 


0.05287 


10.08102 


-0.05814 


0.94518 


0.27414 


7.41921 


0.5 


-0.02668 


0.38076 


0.05430 


8.32835 


-0.05149 


0.78149 


0.28845 


5.93073 


0.6 


-0.02434 


0.32725 


0.06149 


6.75355 


-0.04571 


0.65486 


0.33023 


4.61883 


0.7 


-0.02233 


0.28415 


0.07542 


5.31595 


-0.04026 


0.54723 


0.40290 


3.45588 


0.8 


-0.01862 


0.23208 


0.11519 


3.75473 


-0.03496 


0.44995 


0.50912 


2.45910 


0.9 


-0.02015 


0.22374 


0.11821 


3.07080 


-0.03043 


0.36682 


0.62743 


1.69514 



Note. — bi's are defined through equation H38| l. 



Table 6. Fitting coefficients for Nn of distorted polytropes 







n = 1.5, bo 


= 0.42422 






n = 3, bo 


= 0.36394 




6i 


b2 


b3 


b4 


6i 




bs 


b4 


0.1 


0.00696 


-0.18218 


-0.06256 


14.81447 


0.00110 


-0.06453 


1.14963 


-11.8357 


0.2 


0.00533 


-0.11239 


-0.03379 


11.12954 


0.00086 


-0.03997 


0.56067 


-9.41496 


0.3 


0.00452 


-0.08394 


-0.02569 


8.96846 


0.00075 


-0.03009 


0.36519 


-8.22768 


0.4 


0.00399 


-0.06755 


-0.02301 


7.32661 


0.00068 


-0.02458 


0.26786 


-7.45925 


0.5 


0.00357 


-0.05633 


-0.02307 


5.92820 


0.00064 


-0.02104 


0.20972 


-6.89082 


0.6 


0.00321 


-0.04771 


-0.02524 


4.67990 


0.00061 


-0.01856 


0.17112 


-6.43430 


0.7 


0.00288 


-0.04047 


-0.02957 


3.55301 


0.00059 


-0.01675 


0.14367 


-6.04829 


0.8 


0.00229 


-0.03146 


-0.04086 


2.37704 


0.00058 


-0.01534 


0.12333 


-5.72063 


0.9 


0.00238 


-0.02892 


-0.04344 


1.78659 


0.00055 


-0.01410 


0.10903 


-5.50715 



Note. — bi's are defined through equation I I38I I. 
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Table 7. Fitting coefficients for / of distorted polytropes 



fi 




n = 1.5, bo 


= 90.91016 






n = 3, bo 


= 93.15765 








63 




61 


62 


63 


64 


0.1 


0.23389 


-6.78106 


416.21309 


0.06683 


-4.44876 


133.13261 


152.14256 


10.18132 


0.2 


0.13995 


-3.24888 


197.05415 


0.05320 


-3.27895 


78.57307 


81.88516 


7.48309 


0.3 


0.10835 


-2.25890 


127.38938 


0.00015 


-2.69940 


56.50289 


58.48508 


6.01238 


0.4 


0.09711 


-1.90272 


93.74429 


-0.06256 


-2.32438 


43.95814 


47.13868 


4.97229 


0.5 


0.09662 


-1.80771 


74.20831 


-0.12872 


-2.05058 


35.62331 


40.77503 


4.14597 


0.6 


0.10447 


-1.86636 


61.68649 


-0.19697 


-1.83950 


29.58690 


36.91530 


3.45368 


0.7 


0.12211 


-2.06525 


53.29763 


-0.26863 


-1.68084 


25.05345 


34.30903 


2.86731 


0.8 


0.15798 


-2.48941 


47.93900 


-0.35028 


-1.48905 


20.90588 


33.51635 


2.31784 


0.9 


0.25324 


-3.60727 


46.46436 


-0.47006 


-1.67258 


20.43106 


27.78674 


2.09483 



Note. — bi's are defined through equation I I38I I 



Table 8. Fitting coefficients for fc of the critical components 





n = 1.5 


n = 3 


do 


5.283 


10.051 


di 


22.07 


39.75 


d2 


-87.99 


-155.79 


ds 


236.69 


420.65 


d4 


-292.31 


-522.14 


ds 


141.03 


253.75 


eo 


1.101 


2.242 


ei 


7.922 


14.798 



Note. — di's and e^'s 
are defined through equation 
l l45t . 
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Table 9. Fitting coefficients for AlogA;2,eff 







n=1.5 






n=3 






n=3.5 




ai 


a2 


^3 


ai 


a.2 


13 


ai 




^3 


0.1 


0.73385 


34.10425 


-148.77090 


-0.03012 


24.47968 


-196.59523 


-0.08969 


21.17797 


-182.98519 


0.2 


0.88514 


30.49465 


-120.51368 


0.12028 


21.15117 


-167.29413 


0.07023 


16.73734 


-137.84158 


0.3 


1.00745 


25.75547 


-87.11816 


0.23724 


17.77352 


-135.98160 


0.17016 


14.06672 


-113.16240 


0.4 


1.05899 


23.41619 


-91.19969 


0.32678 


14.69585 


-107.93152 


0.24079 


11.91660 


-93.63182 


0.5 


1.16935 


15.72668 


-29.29898 


0.41656 


10.91861 


-74.44114 


0.29680 


10.04001 


-79.29394 


0.6 
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Note. — ai's are defined through equation JSOj. 
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Fig. 1. — Schematic view of a binary configuration consisting of a distorted polytrope (star I) with 
mass Ml and its point-mass companion (star II) with mass M2. The Cartesian coordinate system 
{x,y,z) corotating with the binary is centered at the center of mass of star I that has polar and 
equatorial radii Rp and Re in the z- and positive x-directions, respectively. Star II is located at 
{A, 0,0). Dotted line draws a sphere with radius equal to the distance from the center of mass of 
star I to the Lagrange point Rl of the binary. 
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Fig. 2. — Dependences on the rotation parameter v of (a) the mass concentration Pc/p, (b) volume 
radius ^i, (c) equatorial radius ^e, and (d) polar radius for rotating polytropes with n = 3. 
Solid lines plot the results of the present s tudy, while dott e d and dot- dashed l i nes ar e for the first- 
and second-order perturbation analyses of IChandrasekhad (Il933al ) and lAnandl ( 196811 . respe ctively. 
Open circles denote the values from the adva nced perturbation theory of lLinnelll (Il977l . ll98ll ). while 
the direct numerical results of iJamesI (|l964l ) are given by filled circles. 
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Fig. 4. — a: Critical values of the aspect ratio TZ = ^e/^p resulting from our self-consistent solutions 
for rotating polytropes with index n {solid line) and the prediction, TZc = 1.5, of the Roche models 
{dotted line), b: Dependence on TZ of the ratio Vg/Vc of the equatorial to critical velocities. Curves 
for n > 3.0 are almost indistinguishable. 



-36- 




Fig. 5. — Dependences of the rotation parameter v on the dimensionless orbital separation v for 
self-consistent equilibrium polytropes with (left) n = 1.5 and {right) n = 3 under both rotational 
and tidal perturbations. Each solid line corresponds to a sequence of equilibrium polytropes with 
fixed relative mass // indicated at one end of the line. The critical configuration at each sequence 
is marked by a filled circle. The dotted line connecting the filled circles demarcates the boundary 
in the v-v plane above which no equilibrium polytrope exists. 
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Fig. 6. — a-c: Meridional cross-sections of constant density surfaces of distorted equilibrium poly- 
tropes with ra = 3. Solid contours correspond to p/p = 3~™'(pc/p)iso with m = 0, 1, • • • , 5, oo from 
inside to outside, d: Outer boundaries of the distorted polytopes in critical configuration. In all 
the panels, the dotted line draws the Roche lobe around star I that has a half of the total mass 
{jj, = 0.5). See text for details. 
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Fig. 7. — Comparison of the volume radius Ri of a critical component with the effective radius Rl 
of the corresponding Roche model. 
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Fig. 8. — Dependences of Pc/p and on the relative radius Ri/A = ^i/u for {left) n = 1.5 and 
{right) n = 3.0 polytropes. Each sohd line is a sequence of distorted polytropes with fixed relative 
mass jjL indicated at one end of the line. The dotted lines connecting the filled circles, denoting the 
end points of the sequences, correspond to the critical components. 
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Fig. 9. — Same as Figure M for Nn and /. 
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Fig. 10. — Changes in the radius of a star with fixed mass due to tidal and rotational perturbations 
under the assumption of constant pressure constant K for {left) n = 1.5 and (right) n = 3.5 
polytropes. Each solid line is a sequence of distorted polytropes with fixed relative mass /x indicated 
at one end of the line. The dotted lines correspond to the critical components. 
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Fig. 11. — Ratio of the spin to orbital angular momenta, Jspin/<^rb) of (left) n = 1.5 and (right) 
n = 3.0 polytropes, as functions of ^i/i' and /v,. Each solid line corresponds to a sequence of distorted 
polytropes with fixed relative mass fi indicated at one end of tlic line. The dotted lines connecting 
the filled circles, denoting the end points of the sequences, represent the critical components. 
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Fig. 12. — Dependences on fi of (a) the critical rotation parameter Vc and (b) the critical dimen- 
sionless o rbital separatio n I'r fo r polytropes with n = 1.5 and 3.0. In (a), the filled circles plots the 
results of ISingh &: Singhl (jl983l ) based on the double approximation method. 
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Fig. 14. — Variations of A log k2 with (left) A and (right) A log p/ pc for polytropes with n = 1.5, 3.0, 
and 3.5. Each solid line is a sequence of distorted polytropes with fixed relative mass /i, with a filled 
circle at the tip of the line c orresponding to a cri tical cornpone nt. For comparison, A log k2 = — 0.7A 
and — 0.87A, the results of IStotherd (jl974l ) and IClaretl (jl999l ) for rotationally-disturbed stars, are 
niotterl a.s dotterl a.nri rla.sherl lines resnertivelv 
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Fig. 15. — Apsidal motion rate w versus the relative radius ^xjv for cases with (lefi) /x = 0.1 
and {right) ji = 1.0. Sohd Unes plot vjn calculated from our self-consistent solutions for distorted 

polytropes, while dotted lines arc for the theoretical rates wt based on the first-order perturbation 
method, with k2 calculated for distorted polytropes from the Radau equation. For sufficiently large 

lij. rn„ is IflrTOr than tttj. 




Fig. 16. — Relative difference between the numerical and theoretical apsidal motion rates, Wn and 
wt, against the relative mass /x for critical components. 
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Fig. 17. — Dependence of AlogA;2,efr on the rotation parameter A for {left) n = 1.5 and {right) 
n = 3.5 polytropes. Each sohd line is a sequence of distorted polytropes with fixed relativ e mass 
/i. Do tted lines correspond to the critical components. The result, Alog/c2cfT = — 0.87A, of IClaret 
]'or rotationally-disturbed stars is drawn as a dashed line in the right panel. 



